Quantum rings for beginners II: Bosons versus fermions 



M. Manninen a , S. Viefers'''*, and S.M. Reimann c 
a Nanoscience Center, Department of Physics, FIN-40014 University of Jyvaskyla, Finland 
b Department of Physics, University of Oslo, P.O. Box 1048 Blindern, 0316 Oslo, Norway 
c Mathematical Physics, LTH, Lund University, P.O. Box 118, S-22100 Lund, Sweden 
'Corresponding author. Email address: s.f.viefers@fys.uio.no (S. Viefers) 
(Dated: October 3, 2012) 

The purpose of this overview article, which can be viewed as a supplement to our previous review 
on quantum rings, [S. Viefers et al, Physica E 21 (2004), 1-35], is to highlight the differences of boson 
and fermion systems in one-dimensional (ID) and quasi-one-dimensional (Q1D) quantum rings. In 
particular this involves comparing their many-body spectra and other properties, in various regimes 
and models, including spinless and spinful particles, finite versus infinite interaction, and continuum 
versus lattice models. Our aim is to present the topic in a comprehensive way, focusing on small 
systems where the many-body problem can be solved exactly. Mapping out the similarities and 
differences between the bosonic and fermionic cases is of renewed interest due to the experimental 
developments in recent years, allowing for more controlled fabrication of both fermionic and bosonic 
quantum rings. 
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I. INTRODUCTION 



Strictly ID systems are unusual in that they provide us with one of the rare cases of exactly solvable models for a 
quantum mechanical many-particle system. Thus, a tremendous amount of work has been dedicated to such systems. 



studied the classical equation of 



We mention a few selected highlights illustrating the history of the field: Tonks 

state of elastical spheres in ID, and the corresponding quantum mechanical derivation was given by Nagamiya 
Girardeauf^, 0] showed that there is an intimate relationship between impenetrable bosons and fermions in ID 



Nowadays, ID gases with infinitely strong short range repulsion thus go by the name Tonks- Girardeau gases. Lieb 



gerjj]] a 

(This topic has later on been taken up by many authors, see for example, Refs. 



the ID fermi gas was solved exactly by Lieb and Wu 
Kolomcisky and Straley 



12 



). The Hubbard model for 



14j for fermions and by Girardeau and Wright 



131 ] . The quantum mechanics in ID has been reviewed by 
for Tonks gases. More recently, Girardeau 



15 



and Minguzzijla] have studied quite generally soluble models for ID gas mixtures, and also discussed the fermionic 
[n]]. Infinitely long ID fermion systems arc usually described as Luttinger liquids [2, [l^. In ID 



Tonks-Girardcau gas [17] . Infinitely long ID fermion systems are usually described as Luttinger liquids 
the Fermi surface consists only of two points, which leads to a so-called Pcierls instability 



2C I and a breakdown of the 



Fermi liquid theory. In such a system the low energy excitations are collective, have a phonon-like linear dispersion 



and thus behave like bosons 



21 



- 231] . Arguably the most exotic property of a Luttinger liquid is spin-charge separation 
- spin and charge excitations of the interacting ID system generally move at different velocities. In fact, in small finite 
quantum rings where the many-body spectrum can be solved exactly, the separation of phonon-like charge excitation 



24, 



25 



and spin excitations, which can be described in the Hcisenbcrg model, is seen in a very explicit manner 

Particles confined in a quasi-onc-dimcnsional ring can support quantum states exhibiting persistent currents in 
normal metals l^. \^/\ . superconductors, quantum liquids and quan tum gases. Persistent currents have mostly been 



studied in normal metals and superconductors (for reviews see 



28 



29(). Superfluidity of helium liquids provides 



strongly interacting atomic systems where persistent currents were predicted 



bosons in 4 Hc 



30 



321 ] and for fermions in 3 He 



33 



30 



31 



34| . The possibility of experimentally realizing annular traps 



and subsequently observed for 

44] 



44 



35 



Such 



for cold atom condensates opened up new possibilities to study persistent currents in bosonic ring traps 
systems have been analyzed in many theory works, most of them applying the Gross-Pitaevskii approximation, but 
also going beyond; see for example, Refs. 



47 



55|. 



In metallic or semiconducting quantum rings the current can be generated by a magnetic flux through a mesoscopic 
ring. For neutral atoms confined in a trap, with electromagnetic fields it is possible to create interactions which 

56l459j. Also, microchip design for persistent current measurements has been 



mimic the effect of a magnetic flux 



suggested 
(see [61 



43 



60j. In addition, there has been a lot of recent interest in atomic quantum gases with dipolar interactions 
621 ] for reviews). In ring traps, there is an interesting interplay between the anisotropy of the dipolar two-body 



interactions and the topology of the ring trap; we do not discuss such dipolar systems any further here, but refer to 
the recent works in Rcfs. 63-671. 



The key difference between bosons and fermions is the symmetry of their many-body wave functions: Exchange of 
the (spatial and spin) coordinates switches the sign of the fermion wave function, but leaves the boson wave function 
unchanged. In general, the spin of bosons is an integer while that of fermions is a half-integer. In the case that 
the Hamiltonian is independent of spin, the spin only adds an additional degree of freedom which has to be taken 
into account in requiring the proper symmetry of the many-body wave function. It is often illustrative to talk about 
pseudo- or isospin for systems with different components of particles. We can have one-component fermions, for 
example, by polarizing an electron gas so that each electron is in the spin-up state. In this case one usually speaks of 
"spinless" fermions, in analogy to bosons with spin zero. In the case of cold atom clouds one can prepare a bosonic 
condensate which has atoms in two different hypcrfine states. This is a two-component boson system, and if we choose 
to denote the different hypcrfine states as different (pseudo) spin states, we effectively have bosons with pscudospin 
1/2. 

The purpose of this review is to highlight the differences of boson and fermion systems in one-dimensional (ID) 
and quasi-one-dimensional (Q1D) quantum rings. (For harmonic traps set rotating, a comparison between fermion 



and boson quantum droplets is given in Refs. 



68 



69( , also discussing the formation of vortices) . 

Mapping out the similarities and differences between bosonic and fermionic rings is of renewed interest due to the 
experimental progress during recent years, in particular concerning trapping techniques for cold atoms. While the orig- 
inal experiments on quantum rings were performed on electronic systems (such as in semiconductor heterostructures), 
bosonic as well as fermionic counterparts can nowadays be fabricated with appropriately trapped (charge-neutral) 
cold atoms. These systems allow for a much larger degree of control of physical parameters (such as interaction 
strength) than is the case for electrons in semiconductor structures. Thus, while a lot of interesting work has been 



70l-l73l| , we here mainly have in mind atomic systems where interactions 



done on Coulomb-interacting electronic rings 
are typically short-range. 

The main emphasis in this paper is put on small systems where the many-body problem can be solved exactly. Our 
intention is to keep the theory as simple as possible, and avoid going into too much detail of advanced many-particle 
descriptions. The present paper provides a natural supplement to our previous review on quantum rings, Ref. 25j |. 
which we will frequently refer back to. 

The outline of the paper is as follows. We first consider one-component boson and fermion systems. We start 
with particles interacting with infinitely strong delta function (contact) interactions. In section |TT] we study the 
exact solution of the two-particle case, and in sections IIIII and IIVI we generalize the solution to many particles and 
show how the boson and fermion spectra are similar for odd numbers of particles, but differ for even numbers of 
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particles. In Section [V] we introduce the concept of persistent currents mentioned previously, and demonstrate the 
relation between the current and angular momentum in small rings. In Section IVTl we discuss lattice models for rings 
with only one component of particles. Section IVIII considers two-component bosons and fermions. In section IVIIII 
we abandon the ideal case of infinitely strong interactions, and study the effect of finite interparticle interactions. 
Quasi-one-dimensional systems are considered in Section IIXI 



II. TWO SPINLESS PARTICLES WITH INFINITELY STRONG CONTACT INTERACTIONS 



We start with the simplest case, namely two spinless particles confined in an infinitely narrow ring and interacting 
with an infinitely strong repulsive contact interaction that has the form of a 5-function, 5(x\ — x 2 ). When speaking 
of "spinless fermions", what one typically has in mind is fully polarized electrons (all in the same spin state, say 
spin-up). The Hamiltonian for such a system is (in atomic units K = m = 1) 



Ho 



i / a 2 



d' 2 



2\dx\ dx\ 



+ 9oo5{x 1 - x 2 ), 



(1) 



where x is the coordinate along the ring. For convenience we choose R = 1 and thus x can be considered as an 
angle between and 27r. We assume the limit — > oo, so that the wave function is zero whenever x\ = x 2 (which 
for fermions is automatically satisfied due to the Pauli principle). With the change of variables u = x\ — x 2 and 
v = (xi + x 2 )/2 we get 



Ho 



( d 2 icP_ 

\ du 2 4 dv 2 



The solution can be written as 



ip(u, v) = B(u) sm(nu/2)e l 



(2) 



(3) 



where B is a piecewise constant function needed for bosons and explained below. The periodic boundary condition 
xf)(u,v) = ^{u^v + 2ir) requires that m is an integer, and the condition ip(0,v) = ip(2ir,v) = requires that n is a 
positive integer (a negative integer would only change the sign of the wave function). Using the original coordinates 
we can write the wave function as 



tP(x u x 2 ) = B(xi - x 2 ) sin gfo - x 2 )) e M«i-H*)/2. 



(4) 



We first consider fermions. In this case the function B is not needed (or B(u) = 1) since the rest of the wave 
function is already antisymmetric with respect to coordinate change. We have an additional requirement of periodicity 
ip[x\, x 2 ) = tp(xi + 2-7T, x 2 ) = ip(x\, x 2 + 27r), which is fulfilled only if n + m is even. It is easy to see that in this case 
the solution is a single Slater determinant 



ip(xi,x 2 ) oc exp 



n + m 
1 — 7. — x i 



exp 



m — n 
i — - — x 2 



exp 



?i — m 
i x\ 



exp 



m + n 
1 — 7. — x 2 



(5) 



consisting of single-particle states with integer angular momenta (m + n)/2 and (to — n)/2. 

The case of bosons is more cumbersome since we have to require that the wave function is symmetric with respect 
to particle exchange. This can be done with a piccewise constant function B which has constant absolute value 
but changes sign at x\ — X2 — 2irk where k is an integer. Note that the discontinuity of B is not a problem since 
at these points the wave function is zero due to the sine function. The required function is then a square wave 
B(u) = sgn(sin(it/2)), where sgn is the sign function. The requirement that ib(x±, x^) = ip{x\ +2ir, x%) is then fulfilled 
for bosons only if n + m is odd, since adding 2?t either to x\ or to X2 changes the sign of B(u). 

Note that the bosonic wave function is not analytic at x\ = X2 (for odd n), but has a cusp. However, this is allowed 
since we assumed that the contact interaction is infinite and, consequently, the wave function has a node at that point. 

In general, we can now write the eigenvalues for the two-particle system with infinitely strong contact interaction 

as 

1 



-(n + to ), 



(6) 



where n + m has to be even for fermions and odd for bosons 
The total angular momentum of the state is to, since 

Lip ■ 9 9 



' dx\ + dx2 ' 



-i—tb = mip. 
ov 



(7) 



The wave function thus consists of an internal part, sin(nu/2), and a center of mass part, exp(imv). The latter 
corresponds to a rigid rotation of the two-particle system, while the internal part corresponds to an internal vibrational 



mode[25(. Note that due to symmetry restrictions, we get the same selection rules as for two-atomic molecules: Purely 
rotational states (n = 1) are allowed for fermions only with odd angular momenta, and for bosons only for even angular 
momenta 



74 1 . Figure [T] shows the energy spectra displaying the alternating allowed fermion and boson states for each 
of the vibrational states. 

The ground state of two bosons (to = 0, n = 1) has zero angular momentum, energy E = 1/4 and, interestingly, it 
naively looks like a function of single-particle states with half-integer angular momenta (or half of the mass), 



x j / 2bosons _ 



sgn 



Xl - x 2 
2 

Xl - x 2 



1 

2i 



Xl - X 2 

2 

ixi/2g — 1x2/2 



-iii/2 J12/2 



(8) 



These, however, are not single-particle eigenstatcs of the ring. Fourier expansion in terms of the true single-particle 
cigenstatcs gives an infinite Fourier scries (denoting, as earlier, x\ — X2 = u) 



,T,2bosons 



2 ^> 4 cos(tou) 
7r ^ 7r(l - 4m 2 ) 

2 

(1 -4to 2 ) 



2 

K 



^ 7T(1 



+ e 



(9) 
(10) 
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FIG. 1. Energy spectra for two spinless fermions (blue) and for two bosons interacting with infinitely strong delta function 
interactions (pink). The lines show the different 'vibrational modes' 

In the case of spinless fermions the antisymmetry dictates that a single Slater determinant is always a solution 
for contact interactions, as pointed out above. The ground state, found by setting n = m = 1 in j5]), equals 
exp(zxi) — exp(ix2)- It has angular momentum 1 and energy E = 1/2, while the state with zero angular momentum 
(to = 0, n = 2) is 

YT>2fcrmions „ixi „ — 1x2 ~ix\Ax2 

w l=o — e e — e e 

and has higher energy E = 1. Both energies are higher than the ground state energy for bosons (E = 1/4), as 
illustrated in Fig. [TJ 

In contrast to spinless fermions, for bosons the interaction has a true effect. This allows the bosonic wave function 
to 'take advantage' of all single-particle states, like in Eq. (fTU|) . for lowering the energy. 



III. SPINLESS FERMIONS WITH CONTACT INTERACTIONS (NONINTERACTING FERMIONS) 

In this section we discuss the so-called yrast states, i.e. lowest energy states for given total angular momentum, of 
non-interacting spinless fermions in a strictly ID ring. Most of the detailed analysis is presented for an odd number 
of particles, with the corresponding results for even N stated at the end. 

The ground state of N noninteracting fermions is a Slater determinant of the lowest TV single-particle states. In a 
strictly one-dimensional (ID) ring with radius R = 1 the single-particle energies and wave functions are 

e< = ~, = e Ux , (11) 

where the I is the single-particle angular momentum and x the coordinate along the ring (equal to the angle in polar 
coordinates). In the ground state the TV lowest states are occupied. For an odd N the single-particle states from £ = 



to £ = ±(N — l)/2 arc occupied and the energy becomes 



E°™(N) = 2 £ L = ±N(N>-1). 



(12) 



The angular momentum of the ground state is zero. The corresponding ground state wave function for an odd number 
N of fermions is the Slater determinant 



^GS(X1,X 2 , ■ ■ ■ ,X N ) 



—i(n — l)xi g—i(n — l)x2 , . , g— i{n— l)xf 
Anx\ AnX2 ~inXN 



where n = (N — l)/2. Using the mathematical properties of a determinant we can recast this as 



* GS = e- WM!l e 



N 



IXi JIX 



where in the last step wc have identified the determinant as a so-called Vandcrmonde determinant 75 1. 

By noticing that each Xi appears in 2n factors of the product (2n = N — 1), we can distribute the exponentials in 
front into these factors and get [3] 



N 



N 



*gs =Y[(e«*<-*M 2 -e***-*^ 2 ) = (2 J ) 7V ^ 1 )/ 2 .[]sin 

i<j i<j 



(13) 



Since we are working with unnormalized wave functions, the overall constant factor will be ignored in the following. 

It is easy to verify that multiplying the wave function of this (or any other) many-particle cigenstate with the 
symmetric function exp(i ^2 x i) will again generate a solution of the many-particle problem, with an energy increase 
E — > E + AI + N/2 where M is the total angular momentum of the original state. This operation corresponds to an 
increase of the angular momentum of each single-particle state by one. Exciting the ground state this way (recalling 
that its total angular momentum is zero) will thus increase the angular momentum by AM = N, and the total energy 
by N/2 = AM 2 /2N = AM 2 /2I. This term in fact has the same form as the classical rotational energy of a rigid 
body, with moment of inertia ./V (we have chosen the radius and mass to be equal to one). For this reason we will 
refer to it, somewhat sloppily, as "classical rotation energy", although of course it results from a purely quantum 
mechanical derivation. More generally, we find that the lowest energy states for total angular momenta M = mN are 



obtained through rigid rotation of the ground state, with wave function 

N N 



iV iV / \ 

4 1V1 \ ^ ' 



and corresponding energy 

1 M 2 
fiW = YA N{N2 - 1)+ 2N> (14) 

where the last term is again recognized as corresponding to rigid rotation of the entire ring. 

For angular momenta different from mN we can also easily determine the lowest energy from the single Slater 

determinant of each state. For example, the lowest state for M = 1 is obtained by exciting the topmost single-particle 

from £ = (N — l)/2 to I = (N + l)/2. In terms of wave functions this means multiplying the ground state with a 

symmetric polynomial of exponentials e~xp(ixj), 



TV \ N 



with energy Em=i = Eqs + N/2. We note that, curiously, the energy is the same for M = 1 and M — N. This energy 

increase, N/2, is much larger than what the classical rotational energy, 1/2N, would be. The energy difference is a 

vibrational energy of the quantum state. 

It is straightforward to calculate the lowest vibrational energy for any angular momentum v (1 < v < M), i.e. the 

energy difference between the true quantum mechanical energy and the 'classical' rotational energy: 

. . v(N + 1 - v) v 2 
AEJv) = — — , (16) 

where the first term is the lowest excitation energy from the M = to the M = v state, obtained by exciting the 
particle at i = (N - l)/2 - (v - 1) to I = (TV - l)/2 + 1 (see Fig.l in Ref. 



251 ) , and the second term is the subtraction 



of the classical rotation energy. Since any state can be multiplied by the rigid rotation exp(z?7i ^2 x{), we obtain that 
the rotational spectrum has 'a period of TV' in the sense that the internal structure of the state does not change when 
the angular momentum is increased by iV 25| . This means that the vibrational energy must be the same for M = u 



and M = v + mN (m integer). The above 'phonon' spectrum, taking the form of an inverted parabola, is analogous 
to a classical row of particles interacting with 1/al 2 interaction 



25 



77 79]. In quantum mechanics the effective 1/d 



interaction results from the kinetic energy: A particle confined between the neighboring particles is essentially a 
particle in an infinite box with kinetic energy proportional to 1/d 2 . Classical particles interacting only with a contact 
interaction would not have vibrational states. 

We can now construct the lowest energy for each angular momentum, i.e. the yrast states, by adding the ground 
state energy for M — 0, i.e. Eq. (|T^1) . the 'classical' rotational energy M 2 /2N , and the vibrational energy of Eq. 
P6p . taking into account that it is the same for v and v + mN . The result can be written as 

1 1 1 M 2 1 

E°»£(N, M) = -{N + 1)(M mod N)--(l+ -)(M mod N) 2 + — + -N(N 2 - 1), (17) 



10 

as may be independently verified by a numerical computation of the quantum mechanical spectrum. Higher 'vibra- 
tional' excited states can be derived from the single-particle energies in a similar manner. 

Finally, let us state the corresponding main results for an even number of spinlcss fermions. In this case the 
ground state can again be expressed as a Slater determinant, but now the angular momentum is not zero but rather 
M = ±N/2. Manipulations similar as those leading to Eq. pU)) give the wave functions for the two degenerate ground 
states 

*§ s mions = n e±Mi/2 n sin { x -^p-] > as) 

where total angular momentum can be immediately read off from the first product. An analysis similar to the one 
leading to Eq. lfTT]) can be carried out to find the yrast energies, 

1 1 1 M 2 1 

E^(N 7 M) = -(N + 1)((M + N/2) mod N) -(1 + -)((M + N/2) mod A) 2 + — + —N(N 2 1). (19) 

IV. BOSONS WITH INFINITELY STRONG CONTACT INTERACTIONS 



For odd numbers of bosons the solutions of energy levels are equal to those of spinlcss fermions. This can be 
most easily shown in a lattice model which we will return to in section IVT1 The ground state wave function can be 
obtained from the fermion solution by multiplying it with an appropriate product of sign functions, to make the state 
symmetric with respect to particle coordinate change. Like in the two-particle case, the boson wave function is not 
analytic, since at the points Xi = Xj the partial derivatives are not continuous. However, this is not required since 
we assume an infinitely strong repulsive interaction (in fact even a finite contact interaction causes a discontinuity 
in the derivative of the wave functionjsj]). For example, the ground state for an even number of bosons is a direct 
generalization of the two-particle case ([5]) , 



JY 



^bosons = JJ sgn 



1<J 

and has angular momentum M = 0. 



Xi Xn 



X <t X/ 



54 . It is qualitatively similar for odd and 



In the boson case the energy of the yrast line was recently studied in Rcf. 
even numbers of particles and is expressed in both cases by Eq. (JTTJ) . This means that for odd numbers of particles 
the yrast lines for bosons and fermions are identical to each other, while for even numbers they differ qualitatively. 
Figure [5] shows the energy spectra for 8 spinlcss fermions and bosons with infinitely strong contact interactions. Like 
in the case of two particles, the boson and fermion states have different 'selection rules' which allow combination of 
the angular momenta of the 'vibrational states'. 

Again we find that the energy spectrum is periodic in TV in the following sense: For any energy eigenvalue E\m\ wc 
have E]M'=M+kN] = E[M] + k 2 N/2, and at the yrast line 2£[m=i] = E1[M=N] = -Egs + N/2. In the case of bosons this 



11 




FIG. 2. Energy spectra for 8 spinless fermions (blue) and for 8 bosons interacting with infinitely strong contact interactions 
(pink) in a ID quantum ring of radius R = 1. The small dots show forbidden states. The blue and pink lines show the yrast 
energies for fermions and bosons. 



holds both for odd and even numbers of particles. 

As already mentioned in the introduction, we here stress again that boson systems with finite contact interaction 
have been solved exa ctly 



example. Refs. 



12 



50 



already long time ago by Lieb and Liniger[5j, |6[ and 7|. and later on many others (see for 



54 1, this list being far from complete). 



V. PERSISTENT CURRENT AND ANGULAR MOMENTUM 



In the case of electrons, or other charged particles, a persistent current can be introduced by a magnetic flux through 
the ring. Neutral particles obviously do not respond to magnetic fields in the same way. However, as mentioned in 
the Introduction, there are ways of mimicking the effect of a magnetic flux on, e.g. charge neutral bosonic atoms 



43 



3- 



56H59j. Thus, the following discussion is relevant also to such systems, although we will stick to the standard 
terminology of 'real' magnetic fields. 

The physics of persistent currents becomes particularly transparent if one chooses the magnetic field at the ring 
radius to zero, while the flux through the ring is nonzero. This can be modeled by an appropriate choice of vector 
potential[25j. In such a system the single-particle wave functions of the ideal ID ring do not change due to the flux, 



while the single-particle energies are modified (for a derivation, see e.g. 

2 



25 



1 



em ($) = -( m -_ 



<I> 



(20) 



where $ is the flux through the ring and $o the flux quantum <f>o — h/e (<f>o = 2-7T in our units where K = e = 1). 
The many-particle state for noninteracting electrons (or spinless electrons with contact interactions) is still the same 



12 

single Slater determinant as without the flux, but the total many-particle energy becomes 

E M ($)=E M (0)-M^ + ^ (21) 

In the case of bosons interacting with an infinitely strong contact interaction, the many-particle state consists of many 
pcrmanents made of the single-particle wave functions. Since the flux does not change the single-particle states and 
each permanent must have the same angular momentum, the flux does not change the many-body wave functions. 
The many-particle energy for bosons will also be described with Eq. (|2ip above. 




contact interact. 
5 10 15 

M 



FIG. 3. Energy spectra for 8 spinless fermions (blue) and for 8 bosons interacting with infinitely strong contact interactions 
(pink) in a ID quantum ring of radius R = 1. A magnetic flux $ = $o passes through the ring. The small dots show forbidden 
states. The blue and pink lines show the yrast energies for fermions and bosons. 

The current (density) for a sin gle-p article state cj>i in a ID system, with a flux $ passing through the ring (of radius 



R = 1 and in units H = e = 1), is 



25] 



* = 5 fal - - + ^mm) = i~w (22) 

It follows that the current of the many-particle state, consisting of Slater determinants or permanents, will be 

J = M/$ — iV$/$o- (23) 

This is because each determinant (permanent) has the same total angular momentum M and consists of N single- 
particle states. 

Generally, it can be shown that the current caused by the flux can be derived from the flux dependence of the total 



energy as 



25 



Clearly, when applied to Eq. (|2Tj) the above equation gives Eq. (j23|) . At finite temperatures the current can be 
derived by replacing the total energy E with the Hclmholtz free energy F in Eq. 



13 




* / $o * / *u 

FIG. 4. Energy spectra for 8 spinless fermions (blue) and for 8 bosons (pink) as a function of the flux passing through the 
ring (left), and the corresponding persistent current (right). Energies and currents corresponding to different angular momenta 
are shown as dashed lines and those corresponding to the ground state as solid lines. The particles interact with an infinitely 
strong contact interaction. 

Comparing figures [2] and [3] we notice that the effect of the flux is to tilt the spectrum so that energy states with 
higher angular momenta become the lowest energy states. When we determine the persistent current for the ground 
state we notice that we jump from one angular momentum state to another when the flux increases. For 8 bosons 
these states are 0, 8, 16 etc., and for fermions 4, 12, 20 etc. (see left panel in Fig. HJ) - The current as a function of the 
flux for 8 particles interacting with infinitely strong contact interaction is shown in the right panel of Fig. [4] 



VI. HUBBARD MODEL WITH INFINITE INTERACTION (i-MODEL) 



A. Tight-binding basis and spinless particles 



The Hubbard model 



81 



25 



82 



has been extensively used for studying persistent currents in quantum rings (see e.g. 
84j). The basis of the Hubbard model is the tight binding (TB) model: The electrons are predominantly localized 
on atoms at lattice sites, but are allowed to hop between neighbouring sites. In a ID lattice the simple TB model 
is, somewhat surprisingly, also closely related to the free electron model. A numerical solution of the free electron 
Hamiltonian in a lattice leads to the simple TB model with one state per lattice site 



This means that the low 



energy spectrum of the TB model is identical with that of the free-electron case. The simple TB Hamiltonian for the 
ring is 



H = -t^2 of tn+i + HC, 



(25) 



where a\ and Oj are creation and annihilation operators for site i, HC denotes Hermitian conjugate, and the periodic 
boundary condition requires that ajr +1 = a\. The corresponding Hamiltonian matrix has nonzero matrix elements 
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Hij = — t only if i and j are nearest neighbours. The eigenvalues for the Hamiltonian matrix for a ring of L sites and 
radius R = 1 arc (assuming L > 2) 

e(n) = -2tcos (^ n ^) i ( 26 ) 

and the corresponding eigenvectors (wave functions at site j) 

Mj) = e l2 ™ l/L - (27) 

From the periodic boundary condition it follows that e(n) = e(L — n) = e(— n) and <j> n {j) = <f>L-n{,~3) = 0-n( — i)- 
We notice that in a large ring, with L ^ n the low energy eigenvalues are those of free electrons, e(n) « — 2£ + f/c 2 , 
with the wave vector /c = 2im/ L. 

We will first consider the many-body problem of spinless particles interacting with an infinitely strong contact 
interaction in quantum rings made of L lattice sites. The infinitely strong repulsive short strange interaction prevents 
two or more particles from occupying the same lattice site. The many-body Hamiltonian describing such a system is 
the Hubbard model 

nn jj 

H = -tJ^a}a j + -J2Mni-l), (28) 

ij ^ i 

where "nn" means that the sum goes only over nearest neighbors, hi = a\ai and U the onsite interaction. Note that 
for spinless fermions the interaction term vanishes since hf = hi- Also in the case of bosons with U — > oo we can 
neglect the interaction term of the Hamiltonian if we restrict the many-particle basis to states with only zero or single 
occupancy in each site. The Hamiltonian then reduces to the tight binding Hamiltonian, Eq. (|25[) . This model is 
sometimes called the i-modcl. 

We describe the many-body states in the localized basis and write the basis vectors as \n± n^n-i •• -til), where 
occupation numbers m are now restricted to be or 1 (also for bosons) and the total number of particles is n, = N. 
The only difference between the bosons and fermions is thus whether commutation or anticommutation relations are 
used for the operators a and . Consider one term of the Hamiltonian, ajaj+x, where i < L , operating on a state 
with n.; + i = 1 and m = 0. Both for fermions and bosons we get 

oJoj+il-.-Ol---) = -4-1 - - - 10- - ->, (29) 

independent of the occupation of the other states. This is because the creation and annihilation operators operate on 
neighbouring positions and, consequently, the phase factors for fermion operators (powers of (-1)) always add to an 
even number. A difference can appear when i = L, since a^a^+i = (V L a\. In this case the operation with a^+i = ai 
will give no phase factor but the subsequent operation with a} L gives a phase (— l)^ -1 for fermions. Clearly this is 
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+1 if N is odd and —1 if TV is even. The periodic boundary conditions thus give 

— 1 • - • 1 ) for even number of fermions 

(30) 

H~ 1 - - - 1 ) otherwise. 

The important result is that if the number of particles is odd, the Hamiltonian matrices for bosons and fermions 
are identical, leading to the same eigenvalues and eigenvectors. However, two things have to be noticed: (i) the 
eigenvectors are the same only if the localized single-particle basis is used for constructing the many-body states 
and (ii) the meaning of the many-body basis vector nin 2 n 3 ■ ■ • n^) is different for bosons and fermions. In the case 
of fermions it is an antisymmetric Slater determinant while in the case of bosons it is a symmetric function called 
permanent. 

Remembering that at low energies the TB model is equal to the free particle model, it is easy now to reconsider 
the particles in a continuous ring studied in sections IIIII and IIVI In the case of odd numbers of particles the energy 
spectra of bosons and fermions are identical, while for even numbers of particles, as shown in Figs. [T]and[2] they are 
still related, but qualitatively different. 




M -4-3-2-1 12 3 4 



FIG. 5. Energy spectrum of three spinless particles in a Hubbard ring with 10 sites and infinite U. The left panel shows the 
lowest energy levels (dots) as a function of the angular momentum M. The yrast levels are connected with a line to guide the 
eye. The spectrum is identical for bosons and fermions. The right panels show the occupation of the TB single-particle states 
of the ring for the ground state (A) and for an excited state (B), blue for fermions and pink for bosons. The fermion states are 
single determinants with only three of the states occupied, while the boson states consist of numerous permanents. 

We could also use the single-particle states of the TB model, Eq. (p?T)) . when constructing the many-particle basis 
for bosons and fermions. In this case all fermion wave functions would be single Slater determinants, as we noted 
above (the contact interaction has no meaning for spinless fermions). However, a typical boson state would be a 
complicated linear combination of several permanents. This difference is illustrated in Figure [5] which shows the 
energy spectrum for three particles in a ring of ten sites (N = 3, L = 10). The spectrum is identical for bosons and 
fermions. The state vectors are also identical in the occupation number representation when the localized basis is 
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used. However, in the TB basis, where the fermion wave functions reduce to single Slater determinants, the boson 
wave functions remain complicated as shown in Fig. [5] where the occupations of the TB single-particle states are 
shown for the ground state and for an excited state. In the case of fermions only three TB states are occupied. In the 
case of bosons all single-particle states have some occupation. Note that the infinitely strong contact interaction leads 
to a significant depletion of the single particle ground state, which in the non-interacting case would be occupied by 
all TV bosons.. 

The rotational spectrum of a Hubbard ring is characterized by a periodicity of L. In the case when TV <C L the 
yrast line also shows a periodicity of TV in a similar way as in the case of continuous rings. This is illustrated in Fig. 
O showing the yrast line for TV = 3 which has minima at angular momenta -3, 0, and 3. Lattice quantum rings thus 
have all the same features as continuous rings when particles with strong contact interaction are considered. 



B. Magnetic flux and current in a lattice 

The effect of a magnetic flux piercing a ring of L lattice sites is to add a phase in the hopping term of the TB 



Haniiltonian 
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L 



In the lattice we can define the current operator between two points i and j generally as 

2tt 



(/..,<"" «>, - /,.<'"' «>,) ' ( 32 ) 



where is the hopping parameter and <!>,., the phase change between the two points. In a simple ring $^+1 = 
27r$/($o^) = an d the current is the same between any neighbouring sites. Alternatively, we can again 

determine the current as the derivative of the total energy, using Eq. (|24[) . Note that when the ring is quasi-one- 
dimcnsional, Eq. (|32j) is still valid, but for calculating the total current one has to sum the currents through different 
paths. 



C. Particle-hole symmetry 

In a Hubbard ring the persistent current can be determined as the derivative of the total energy as a function of 
the flux (Eq. (|24[) or by computing the expectation value of the current operator between two adjacent points, Eq. 

In the case of bosons with infinitely strong contact interaction (U — > 00) the Hamiltonian is the same for particles 
and holes, i.e. H = ^ajcfj = X] a i a ! since the operators commute when i 7^ j. This means that the ground state 
energy and the persistent current are symmetric with respect to particles and holes, irrespective of the symmetry of 
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the single-particle spectrum. The situation is different for fermions due to the anticommutation rule, which changes 
the sign of the Hamiltonian for holes. Consequently, in the case of fermions the many-body energy and the persistent 
current of the ring are symmetric with respect to particles and holes only if the single-particle spectrum is symmetric. 

VII. TWO-COMPONENT SYSTEMS (SPIN-1/2) 
A. Hubbard model with infinite U 

Let us now now consider particles with spin (e.g. normal spin for fermions and two different hyperfine states for 
bosons), but still apply the same simple Hubbard model with U —> oo. The spin increases the size of the Fock space 
since all spin configurations should be included. Since we still assume the infinite U-, or t-model, only one particle is 
allowed in a lattice site, independent of its spin. The energy spectrum for bosons and fermions is the same if both of 
the spin components have an odd number of particles, or if the total number of particles is odd. In the former case we 
can arrange the single-particle states as |nif, U2f, ■ ■ ■ , n>Lf, nij,, ■ ■ ■ TiLf)- Acting with the fermion operator aj a ai+i l(7 
does not change the sign. In the latter case we arrange the basis states as |nif, 7114., n2f, ri2f, ■ ■ . , n^f, n>Lf), and again 
a\ a ai+i. a can not change the sign (note that only one of rii+i t f and can be nonzero since only one particle is 

allowed in any site. 

The case when both Nf and Nf are even, is more complicated. However, also in this case all the eigenvalues are 
still the same for bosons and fermions, but the degeneracies of the eigenvalues can differ. 

In general, for a fixed total number of particles, N — Nf + JVj,, all the energies of the spectrum are the same for 
any combination of nonzero Nf and Nf, but the degeneracies can differ. To prove this, we first notice that a state 
with Nf = N± (or Nf = N± + 1 for odd number of particles) has S z = (or S z = 1/2), where S z = Nf - Nf is the 
^-component of the total spin. Since the Hamiltonian is independent of the spin, the spectrum of S z = includes all 
energies for states with S z = —N, — N + 1, • ■ • , N, i.e. for all possible values of S z = Nf — Nf. We then only have 
to show that Nf — N — 1 and Nf = 1 already gives all possible eigenvalues. This can be done with the help of the 
Bcthe ansatz solution, as will be discussed in the case of an infinitely strong contact interaction in the section below. 

It is interesting to observe, that adding to a one-component system just one particle of the other component will 
change the energy spectrum completely by allowing now all the states, which were forbidden by the wave function 
symmetry in the one-component case. For example, in the eight particle case of Fig. [5] all the energies shown with 
black dots would become allowed for both bosons and fermions if the 8 particles consisted of, say, 7 " spin- up" particles 
and one "spin-down" particle. 
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B. Bethe ansatz solution for particles with infinitely strong contact interactions 

The Hubbard model for fcrmions in a ID lattice was solved exactly by Lieb and Wu[l^] using the Bcthc ansatz. 
The total energy of a given many-body state with N particles can be written as 

N 

E = -2t^cosfc J , (33) 



where the numbers kj are found by solving a set of algebraic (transcendental) Bethe equations (see e.g. 
limit of infinite U these numbers will be 



25|). In the 



k 3 L 



J N 



(34) 



with 

AT* 



P 



J2 J «- ( 35 ) 



The quantum numbers Ij and J a are related to the charge and spin degrees of freedom, respectively Their possible 
values depend on whether N and Nf are even or odd. Here we will only consider even N, in which case the one- 
component (iV-f = or Nf = N) boson and fermion systems have different spectra. In this case the Ij-.s arc integers 
(half-odd integers), and the J Q :s are half-odd integers (integers) if is even (odd). We note that p will always be an 
integer. Thus, in practice, all eigenvalues for < N-f < N can be found by letting p run over all integers and choosing 
all possible sets of Ij:s with the restriction J max — / m i n < L. 

The total angular momentum depends on the above quantum numbers and can be determined as 

N N t 

M = E^ + E J - ( 36 ) 

3=1 «=1 

For Nf = the quantum numbers Ij are integers and p = 0. In this case we see immediately that the many-body 
energy (f3"3"]) is a sum of the single-particle energies ([25)1 . When L — > oo we recover the result for particles in a 
continuous ring and, say for 8 particles, the spectrum shown with large blue dots in Fig. [2~] Next, consider the case 
JVf = 1. The number p already goes through all integers and, as predicted in the previous subsection, we thus get all 
possible energies, shown in small blue dots in Fig. [21 Increasing jVf from one thus can not give any new energy levels. 
A similar analysis for odd N also shows that with JVf = 1 we already get all possible energy levels. 

We have here considered fermions. In the case of bosons the Hubbard model has turned out to be more difficult 
to solve with the Bethe ansatz method 87 1 . Fortunately, when we are only interested in the energy spectrum of a ID 
ring with infinitely strong contact interaction, we do not need to solve the boson problem, since the fermion solution 
gives the same energies as long as < Nf < N, as discussed in the previous section. 
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C. Noninteracting particles 

The single-particle spectra of noninteracting particles were given by Eqs. (jllj) and ([26]) for a ID continuous ring 
and a TB lattice ring, respectively. In the case of bosons the number of components (number of spin components) 
does not change the energy spectrum, since each single-particle state may be occupied by any number of bosons even 
in the single component case; only the degeneracies of the energy levels depend on the number of components. At zero 
angular momentum all the bosons occupy the I = single-particle state. When the angular momentum is increased, 
the angular momentum state 1=1 starts to be occupied until all particles are in that single-particle state. Then we 
start to occupy the £ = 2 state and so on. The yrast energy increases linearly with the total angular momentum M, 
but the slope has discontinuities at points M = nN , where the next single-particle angular momentum starts to be 
occupied. Figure [6] shows the beginning of the yrast spectrum for eight noninteracting bosons (as pink dots). 




FIG. 6. Left: Yrast spectra of eight noninteracting particles in a two-component system. The lowest (pink) curve shows 
the results for bosons: It only depends on the total number of particles (Nf + N±). The black curve shows the result for 
one-component fermions. The upper red curve represents Nf = 1, the green curve corresponds to Nf = 2, the blue curve to 
Nf — 3, and the lower red curve (small dots) to Nf = 4 (note that the blue and red curve coincide except at two points). 
The small black dots show the result for bosons and fermions with infinitely strong contact interaction for < Nf < 8. Right: 
Dependence of the yrast spectrum of noninteracting fermions on the number of components (spin degrees of freedom). The 
one-component system is shown with big black dots, the two-component system with red dots, the four-component system with 
small black dots, and the eight-component system with pink dots. Note that in the case of bosons the result (pink dots) is 
independent of the number of the components. 

In the case of fermions the number of components has a marked effect on the spectrum of noninteracting particles. 
It is easy to construct the spectrum using the Pauli exclusion principle: Each single-particle state can accommodate 
only one (or zero) particle of each component. Figure [5] shows the yrast lines for 8 particles. The left panel shows 
the result for a two-component system (spin- 1/2 fermions) with different numbers of spin-up and spin-down particles. 
The lowest energies are obtained when the number of spin-up and spin-down particles is the same. The left panel 
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shows the effect of the number of components on the yrast spectrum of fcrmions. When the number of components 
(spin of the particles) increases, the energy approaches that of bosons. When the number of components equals that 
of the number of particles, the fermion spectrum is identical to that of the bosons, since in this case N fermions can 
occupy the same single-particle state. 



VIII. INTERACTION WITH FINITE STRENGTH 



So far we have only considered the ideal limits of zero and infinite (contact) interaction. The more general case 
where the interaction strength is finite is naturally more complicated. Nevertheless, also in this case several exact 
solutions exist. The seminal paper for bosons was that of Lieb and Liniger 5j who solved exactly the problem of 
interacting bosons in ID, and for fermions that of Lieb and Wuflij] who solved the ID Hubbard model for spin- 
1/2 fermions. Instead of reviewing this important analytic work, we here choose to merely illustrate the qualitative 
differences between bosons and fermions, based on exact (numerical) solutions of small Hubbard systems. 



A. Bosons with finite interaction 



In previous sections we have encountered examples of energy spectra for quantum rings of bosons with both infinite 
(Figf2|) and zero interaction strength (FigjB]). The intermediate problem of bosons interacting with finite contact 
interaction was solved exactly by Lieb and Liniger[5|. Naturally, when the contact interaction increases from zero, 
the corresponding yrast lines smoothly interpolate between those of zero and infinitely strong interaction. 




FIG. 7. Dependence of the yrast spectrum on the interaction strength U of the Hubbard model for one-component bosons. 
The ring has 8 sites and 4 bosons. The different curves represent different strengths of the interaction U. From bottom to top 
U — 0, 0.3, 1, 3, 7, 20, and U — > oo. The results for U — and U — !> oo are joined with solid lines, the others with dashed lines. 



To illustrate this, consider a small Hubbard ring with eight sites (L = 8) and four particles (N = 4). Figure [7] 
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shows the yrast energy as a function of the angular momentum for different values of the interparticle interaction U 
of the Hubbard model. The energy is a periodic function of the angular momentum M due to the finite number of 
sites. Nevertheless, even in this small (lattice) ring we already see clearly the same features as for a smooth ID ring: 
In the case of nonintcracting bosons the energy increases linearly with the angular momentum, and in the case of 
U — ^ co the energy has minima at M = and at M = N. When U increases from zero, the yrast energy smoothly 
approaches that for U — > co. The local minimum at M = 4 is reached with rather small interaction strength U = 3. 

In the case of one-component (spinless) fermions the strength of the contact interaction is irrelevant, since the Pauli 
exclusion principle prevents the particles from "seeing" the interaction. 

B. Yrast-line of two-component bosons and fermions 

Two-component systems have an effective spin 1/2. In the case of fermions it can be the true spin of the particles, 
while in the case of bosons it can be the pscudospin, realized for example by two different hyperfine states of bosonic 
atoms. While a contact interaction has no effect for spinless fermions, it does play a role when both spin components 
are present, since the Pauli principle does not prevent particles with opposite spin from occupying the same site. We 
learned above that in the case of an infinite interaction the energy spectrum is the same for bosons and fermions as 
soon as we have at least one particle of the minority component. This not true for interactions of finite strength. 




M M 

FIG. 8. Dependence of the yrast spectrum on the interaction strength U of the Hubbard model. The ring has 8 sites and 4 
particles, 2 spin-up and 2 spin-down. The left panel shows the results for fermions and the right panel for bosons. The black 
points connected with solid lines show the results for the noninteracting case (U = 0), the blue points connected with solid 
lines correspond to U — > oo. Note that this result is the same for bosons and fermions for reasons discussed in Section [VII Al 
The upper curves with dashed lines show results for one-component systems. The boson result for the noninteracting case is 
plotted also in the left panel with pink points connected with dashed lines. The other plots in the fermion case correspond to 
U = 1 (red), U = 3 (blue), U = 7 (green), and in the boson case U = 0.3 (pink), U = 1 (red), (7 = 3 (blue). 

We will first consider the yrast spectrum for a Hubbard system with eight sites and four particles so that we 
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have two particles of each component. Figure [8] shows the results for several strengths of the interaction U. In 
the case of U — > oo the result is the same for fermions and bosons for reasons described in Section IVH Al In the 
case of nonintcracting systems the fermion and boson results differ, except at M — 2 which is the ground state for 
nonintcracting fermions. When the interaction strength increases the yrast lines smoothly approach that of infinite 
U . In the case of fermions the local minimum at M = 2 remains to quite high values of U, it is still marked at U = 7. 
Only at very large U the boson and fermion results will be qualitatively similar. 

C. Spin-splitting of the yrast line 

In the extreme case of infinite U, the effective Hubbard Hamiltonian is totally spin-independent, and thus the 
spectra are degenerate in the spin degree of freedom. As we shall see in this section, this is no longer the case for 
finite interaction strength U, where this degeneracy is split. First, however, let us return to the spin-independent 
limit for a two-component system, i.e. particles of (pseudo)spin 1/2. Different choices of numbers of particles in the 
up- and down components thus correspond to different ^-components of the spin. For each total spin there will be 
a multiplet of states, with the same energy and orbital angular momentum, corresponding to different z-componcnts 
of the total spin or, equivalcntly, different numbers of particles in each component. In the case of four particles, for 
example, the energy spectrum computed for two spin-up and two spin-down particles (iV-j- = 2, ATj, = 2) also includes 
all energy levels for Nf = 3, and iVj, = 1 as well as for N^- = 4 and N± = 0. 
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FIG. 9. Spin-dependence of the yrast spectra for four spin-1/2 fermions (solid lines) and (pseudo)spin-l/2 bosons (dashed 
lines). The results are given for the Hubbard model with 8 sites and U = 20. Different colors indicate the total (pseudo)spin: 
Red for 5 = 0, blue for S = 1 and green for S — 2. The dashed and solid red lines coincide. For infinite U (dashed line with 
dots) the spin-splitting vanishes. For fermions the splitting at large but finite U corresponds to that of the antiferromagnetic 
Heisenberg model. 

The distribution (iV^, N±) also provides information on the total spin of a given state of the two-component system. 



N=4 (pseudospin 1/2) 



■ fermions 
> bosons 
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In the case of four particles we know that the total spin can be 5 = 0, 1, or 2. The totally polarized case with S z = 2, 
i.e. the one-component case, corresponds to 5 = 2. Note that the case = 3, and iVj, = 1 has the z-component of 
the spin S z = 1, but the total spin can be 1 or 2. Similarly, in the case of Nf = 2 and N± = 2, 5 2 = 0, while the total 
spin can be 0, 1 or 2. 

In the case of an infinitely strong interaction, states with different total spin will be degenerate. Away from this 
limit, however, this is no longer true. In particular, in the case of fermions the Hubbard model approaches the 
antifcrromagnetic Heisenberg model for large but finite U 88| , so the Hamiltonian will contain a spin-dependent term 



J N 

Hj = -J2s i -S j . (37) 



2 



. Thus, while 



Here, the effective Heisenberg coupling is proportional to the inverse of U and vanishes at U — > oo 
the energies are independent of the total spin of the quantum state in the infinite U limit, this degeneracy will split 
for finite U. Note that this splitting is caused by the (spin- independent) inter-particle interaction U. 

Figure [9] shows the spin splitting of the lowest energy states for four particles in a Hubbard ring with eight sites. 
The dashed line with dots shows the result for U — > oo. This result is the same for bosons and fermions and naturally 
does not show any spin splitting. The colored lines show the energy levels for different spin states in the case when 
U — 20. In the case of fermions (solid lines) the spin splitting corresponds to that of the antifcrromagnetic Heisenberg 
model. The lowest state for M = has 5=1. This is a result of Hund's first rule: The lowest single-particle state 
(£ = 0) is filled with two particles with opposite spins, while the two remaining particles occupy the states £ = 1 
and £ = — 1 and prefer to have the same spin for removing the effect of the repulsive interparticle interaction. The 
energy of the 5 = 2 state for fermions equals that for infinitely strong interaction. This is because now the system is 
polarized and thus a one-component system where the interaction does not have any meaning. 

The result for bosons is qualitatively different. Now the lowest state has 5 = 2, corresponding to a one-component 
system. The results for 5 = are the same as for fermions. This is true for any value of U as seen also in Fig. [5] 
where all the points for M = 2 are the same for bosons and fermions. In the case of four particles this result appears 
to be independent of the length of the ring. However, computations for six particles show that this is not a general 
property of bosons and fermions in Hubbard rings. 

A major difference between the fermionic and bosonic cases lies in their respective magnetic states. It has been 
shown that the Bose-Hubbard model for spin-1/2 particles with large on-site repulsion can be approximated by a 
ferromagnetic Heisenberg model ] 



90l ] . Again for half-filling the coupling constant has the same value as in the case of 
fermions (J = —4t 2 /U) but now the sign is opposite. In the case of four particles in a ring, the Heisenberg model 
(independent of the sign of J) is exactly solvable (see an exercise in[60]). We notice that the splitting seen in Fig. 
9 is already qualitatively the same as in the Heisenberg model, and that the boson and fermion cases have different 
magnetic order. Numerical diagonalization of the Hubbard and Heisenberg Hamiltonians for four and six fermions 
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and bosons shows that the coupling constant for large U is the same for the antifcrromagnctic fcrmion case and the 
ferromagnetic boson case for all filling fractions [N < L). For half-filling (N/L = 1), \ J\ « At 2 /U , while for smaller 
fillings | J | decreases with N/L[4], 



D. Fixed points in fermion spectra 

Considering particles with infinitely strong contact interactions we learned that in a two-component system it is 
sufficient to have only one particle of the minority component in order to obtain the entire rotational spectrum for 
any mixtures, while the degeneracies will depend on the numbers of particles in each component. In the case of one- 
component (spinlcss) fcrmions the contact interaction does not have any effect and is equivalent to the noninteracting 
case. At angular momenta M = N/2 + kN for even number of fermions and at M — kN for odd number of 
fermions (k any integer), the lowest energy is the same for spinless fermions and for two-component systems with 
N\. = 1 interacting with contact interactions. Moreover, it is easy to discover that the same energy results also for 
noninteracting two-component fermions with iVj, = 1, as seen in Fig. [5] 
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FIG. 10. Dependence of the yrast spectrum on the interaction strength U of the Hubbard model in the case with N = 4, 
N-f = 3, N_i = 1 and L = 8. The left panel shows the results for fermions, and the right panel for bosons. The curves from 
down to up correspond to U — 0, 1, 3, 7, and U — >• oo, respectively. Note that in the case of fermions the energy at M = 2 and 
M — 6 is independent of U. 



This gives us the interesting result that at these angular momenta (M = N/2 + kN or M — kN), as long as N± = 1, 
the lowest energy is completely independent of the strength of the contact interaction. Naturally, this result only 
holds for fcrmions as demonstrated in Fig. 1101 showing the results for N — 4, L = 8 for bosons and fcrmions (now 
7V t = 3 and N± = 1). 
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IX. QUASI-ONE-DIMENSIONAL RINGS 

In experiments, as they were mentioned in the Introduction, the quantum rings will typically not be strictly 
one-dimensional, which means that the particles eventually may "pass" each other. Considering the ring as a quasi- 
one-dimensional (Q1D) wire, the perpendicular modes of the single-particle wave function can be separated from the 
longitudinal modes. If the excitation energy of the perpendicular modes is large, the many-particle state is composed 
mainly of the lowest perpendicular mode, and the system becomes rather similar to the strictly ID case. 



A. Continuous Q1D rings 



A generic continuous Q1D ring can be modelled using a harmonic confinement, 



V(r, z) = -mw 2 (r - R) 2 + V z 



(38) 



where R is the radius of the ring, r the radial coordinate on the plane of the ring, u> the planar confinement frequency 
and m the mass of the particles. V z indicates the confining potential perpendicular to the plane. Depending on 
the system considered, it might be harmonic or, e.g., correspond to the effective confinement experienced by the 2D 
electron gas in semiconductor hetero junctions. Frequently V z is assumed to be so large that motion in the z-direction 
is completely frozen out, and thus the system is effectively two-dimensional. We will here only consider 2D systems 
and thus the second term in Eq. (|38|) is not needed. 

Exact diagonalization techniques as well as density functional methods have been used to describe electrons in 



Q1D rings 
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92(. For reviews we refer to [25 



93j. Exact diagonalization studies show that in narrow rings the 



spin and charge excitations separate: The former resemble those of an antiferromagnetic Heisenberg model, and the 
latter vibrational modes of localized electrons l^. The results of density functional methods (making use of the 
local-spin-density approximation) show clearly the electron localization in narrow rings and give qualitatively correct 



persistent currents 
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97] 



Bosons in annular traps have been studied using the Gross-Pitaevskii m ethod (see for example 



as well as with exact diagonalization (see for example. 
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Figure [TT] shows the yrast spectrum for 17 bosons in a Q1D quantum ring, calculated by Bargi et al. 
diagonalization techniques. In the case of a one-component system (red line) the result is qualitatively similar to that 
of particles interacting with an infinitely strong contact interaction (see Fig. In the two-component case, with 
iVf = 15 and JVj, = 2, a smoother curve is produced. In an infinitely narrow ring the black curve would be expected 
to be a parabola, for the reasons explained in Section IVH Al independent of the number of particles in the minority 



component N± (as long as it is nonzero). It is interesting to note that the effect of the finite width of the confinement 
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FIG. 11. The yrast energy as a function of the angular momentum for a Q1D ring with 17 bosons. The red line is for a 
one-component system and the black curve for a two-component system with 15 particles in the majority component and 2 
particles in the minority component. The energy is in units of the effective oscillator strength of the confinement. From 9ll ]. 

appears to be (qualitatively) stronger in the case of a two-component system than in the case of a one-component 
system. 



B. Q1D rings in the Hubbard model 

Quasi-one-dimensional quantum rings can also be made out of optical lattices which can be described with the 
Hubbard model. In the case of bosons the simple ID Hubbard ring with finite U already mimics the Q1D ring, since 
a finite U means that the particles can pass each other like in a Q1D continuous ring. This is illustrated in Fig. 1121 
The ID ring with a finite U shown in the middle panel of the figure, gives a qualitatively similar low energy spectrum 
as a Q1D ring with infinite U, shown in the lowest panel. 

In the case of fermions the situation is different . In a strictly ID one-component system the Pauli exclusion principle 
prevents particles from passing each other even with zero U . However, in the case of two-component fermions the 
finite U Hubbard model also mimics a Q1D ring '. 



25j . This is illustrated in Fig. [T3] which shows results for seven 
two-component fermions in a ring with 15 sites with the interaction strength [7 = 5. We know from the previous 
sections that for U — > oo the yrast line is the same for any < Nf < N. Figure Q2] shows that this is not the case 
for finite U, but results for different Nf differ qualitatively. It is interesting to note that the case Nf = N — 2 is 
qualitatively similar to the related case for bosons in a Q1D ring, shown in Fig. 1111 
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FIG. 12. Flux dependence of the lowest energy levels for four bosons in quantum rings. The uppermost panel shows the result 
for U — > oo. The middle panel shows the result for U = 3.4. the lowest panel shows the result for U — > oo in the case of a Q1D 
Hubbard ring. The rings in each case are shown at the left. 

The finite width of the quantum ring can change the energy spectrum markedly. In the case of four fermions in a 
ID ring the ground state is degenerate and has a finite angular momentum. This is not necessarily the case if the 
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FIG. 13. Energy as a function of angular momentum for two-component fermions in a Hubbard ring with L = 15 sites, 
Nf + N_i = 7 particles and the interaction strength U = 5. Nf = 7 black; Nf = 6 green; Nf — 5 blue; iV-f = 4 red. For 
comparison, the dashed line shows the U — > oo result which the same for < Nf < 7. Note the qualitative similarity of blue 
curve (Nf — N — 2 = 5) with the black curve in Fig. 1111 



ring is wide enough. The upper panel of Figure [14] shows that in a Q1D Hubbard ring the lowest energy state of the 
fcrmion spectrum is not degenerate and the energy increases with the flux, even if the system has four electrons. In 
the ID case the ground state would be degenerate at zero flux, and increasing the flux would split the degeneracy, 
decreasing one of the energy levels. 

Figure [14] also shows the result for a so-called Mobius ring, studied in detail by Ballon et al. 



101 



where the fcrmion 



ground state is degenerate. It is interesting to note that in the Mobius ring the boson spectrum is nearly identical 
with that of the normal Q1D ring, while the fermion spectrum becomes qualitatively different. Note that only the 
periodic boundary condition is changed in going from the normal ring (upper panel in the figure) to the Mobius ring 
(lower panel). 



X. CONCLUSIONS 



The theoretical description of quantum rings is a vast field, involving various types and strengths of interactions, 
continuum- versus lattice models, strictly ID versus quasi-lD systems, presence or absence of magnetic fields (and 
thus persistent currents), spinful versus spinless particles - and of course the quantum statistics of the particles. Our 
hope is that the present paper has provided a comprehensive and fairly self-consistent overview of the similarities 
and differences between bosonic and fcrmionic quantum rings in all of the above contexts. As we have seen, in many 
cases the quantum statistics of the particles docs make a big difference to the energy spectra, and thus the physical 
properties, of the system. Since both bosonic and fcrmionic quantum rings can be fabricated in the lab nowadays, 
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FIG. 14. Flux dependence of the lowest energy levels for one-component fermions (red) and bosons (blue) interacting with 
infinitely strong contact interactions for the different Hubbard rings shown left. In each case the system has four particles. The 
lower case is the Mobius ring. 



these differences are important to be aware of. The list of references certainly does not exhaust all of the work done 
in the field, in particular since we have chosen to mainly focus on small systems. 
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